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At finite hyperconductivity and finite fermionic density the flux lines of long range hypermag- 
netic fields may not have a topologically trivial structure. The combined evolution of the chemical 
potentials and of pseudoscalar fields (like the axial Higgs) , possibly present for temperatures in the 
TeV range, can twist the hypercharge flux lines, producing, ultimately, hypermagnetic knots (HK). 
The dynamical features of the HK depend upon the various particle physics parameters of the model 
(pseudoscalar masses and couplings, strength of the electroweak phase transition, hyperconductivity 
of the plasma) and upon the magnitude of the primordial flux sitting in topologically trivial con- 
figurations of the hypermagnetic field. We study different cosmological scenarios where HK can be 
generated. We argue that the fermionic number sitting in HK can be released producing a seed for 
the Baryon Asymmetry of the Universe (BAU) provided the typical scale of the knot is larger than 
the diffusivity length scale. We derive constraints on the primordial hypermagnetic flux required 
by our mechanism and we provide a measure of the parity breaking by connecting the degree of 
knottedness of the flux lines to the BAU. We rule out the ordinary axion as a possible candidate 
for production (around temperatures of the order of the GeV) of magnetic knots since the produced 
electromagnetic helicity is negligible (for cosmological standard) if the initial amplitude of the axion 
oscillations is of the order of the Peccei-Quinn breaking scale. 



I. FORMULATION OF THE PROBLEM 



^^ , In a globally neutral plasma, the transversality of the magnetic fields implies, at high conductivity, the conservation 

^s1 ' of both magnetic flux and magnetic helicity. The plasma element evolves always glued together with magnetic flux 
\^ ' lines whose number of associated knots and twists does not change . The energetical and the topological properties 
^^ I of the magnetic flux lines are conserved in the superconducting (or ideal) regime. Moreover, the linearity of the 
O^ ' magnetohydrodynamical equations in the mean electric and magnetic fields implies that there are no chances that 
large scale magnetic fields could be generated only using electromagnetic plasma effects. 

As noticed long ago KM magnetohydrodynamics (MHD) represents a quite powerful method in order to analyze the 

QhI dynamics of the interstellar plasma. In particular the topological properties of the magnetic field lines and of the bulk 

f-^ ■ velocity field play a significant role in the dynamo theory (often invoked in order to explain the typical amplitudes 

q_) ^ of large scale magnetic fields in spiral galaxies 0) and it seems well established that for successfully implementing 

rS^ • the dynamo mechanism we are led to postulate a global parity violation at the scale at which the dynamo action is 

►jj ' turned on |jj§]. Again, because of the linearity of MHD, in order to explain the galactic magnetic fields ||^,|^ we have 

• ^H , to assume that at some finite moment of time (after the Big-Bang) some magnetic seeds M had to be present. 

/\ ' Up to now, the attention has been mainly focused on the energetical properties of the seeds and very little has been 

'^ . said concerning their topological features. Apart from the dynamo mechanism, the topology of the magnetic flux lines 

■ - - ' seems less relevant (for the seed problem) than their associated energy density. However, going backward in time, we 

find that the topological properties of long range (Abelian) gauge fields can have unexpected physical implications. 

The importance of the topological properties of long range (Abelian) hypercharge magnetic fields has been indeed 
emphasized in the past MM . In M the interesting possibility of macroscopic parity breaking at finite fermionic density 
and finite temperature was studied and it was suggested that if parity is globally broken an electromagnetic current 
density directed along the magnetic field can arise. A similar observation has been made in the context of the standard 
electroweak theory at finite temperature and chemical potential |g] where it has been observed that, thanks to the 
anomalous coupling of the hypercharge the fermionic density can be converted into infra-red modes of the hypercharge 
field. 

If the spectrum of hypermagnetic fields is dominated by parity non-invariant configurations the baryon asymmetry 
of the Universe could be the result of the decay of these condensates ||^. The purpose of the present paper is to 
connect the topological properties of the hypermagnetic flux tubes to the Baryon Asymmetry of the Universe (BAU) 
by discussing a class of dynamical mechanisms leading to the generation of HK. 
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The evolution of mean hypermagnetic fields in an electroweak plasma at finite conductivity and finite fermionic 
densities can be derived ||lO[ and it turns out that a good gauge-invariant measure of global parity invariance is 
the hypermagnetic helicity, i.e. Tiy ■ V x Tiy- If the hypermagnetic field distribution is topologically trivial (i.e. 
{Tiy ■ V X Hy) = 0) matter-antimatter fluctuations can be generated at the electroweak scale |ll|. Provided the 
typical correlation scale of the magnetic helicity is larger than (or of the order of) the neutron diffusion scale blue- 
shifted at the electroweak epoch [i.e. L„(re^) ~ 0.3 cm at Tew ~ 100 GeV], matter-antimatter domains can survive 
until the onset of Big-Bang nucleosynthesis (BEN) leading to non standard initial conditions for the light element 
abundances calculations. In the presence of a fraction of antimatter at BBN a reduction of the neutron to proton 
ratio can be foreseen with consequent reduction of the ''He abundance. This effect was recently pointed out in BBN 
calculations ||l2|| in the presence of spherical antimatter domains and further work is in progress |l3| ] . 

An example of topologically trivial distribution of hypermagnetic fields leading to some of the above mentioned 
effects is a stochastic bacground whose (parity-invariant) two-point function can be expressed, in Fourier space, as 

G.,{k)^k'f{k)6^'Hk^k')I^S,,~^y m = ^ (1.1) 

{A{k) is a dimension- less factor measuring the intensity of the stochastic hypermagnetic background). Since hyper- 
magnetic fields are strictly divergence free (i.e. no hypermagnetic monopoles are present) the magnetic flux lines form 
a mixture of closed loops evolving glued to their electroweak plasma element traveling with velocity v. As far as we 
are concerned with scales larger than the hypermagnetic diffusivity scale (i.e. L^^ ^ 3 x lO^^cm at Tew), then these 
flux loops cannot break or twist. 

If the Universe is flUed, at some scale, by topologically non-trivial HK the flux lines form a network of stable knots 
and twists which cannot be broken provided the conductivity is high. Therefore (TYy • V x Tiy) 7^ 0. Thus, completely 
homogeneous hypermagnetic conflgurations cannot be topologically non-trivial. If we want to have knots and twists in 
the flux lines, we are necessarily driven towards field configurations which are not only time but also space dependent 
like the Chern-Simons waves [||£ 
analysis of dynamo instabilities 1 5 



analogous to the magnetic knot conflgurations [[10[ frequently employed in the 
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If the hypermagnetic configurations are topologically non-trivial (and consequently inhomogeneous) some knots in 
the lines of the velocity field can be also expected. Indeed, by solving the hypermagnetic diffusivity equation to 
leading order in the resistivity expansion, the velocity field turns out to be proportional to the hypermagnetic field, 
namely vy — Tiy / {^J N^f fT^^) ■ Notice that this is nothing but the electroweak analogous of the Alfven velocity llq] 
so important in the spectrum of (ordinary) plasma waves |16| . Since the hypermagnetic fields are divergence- free, the 
fluid will also be incompressible (i.e. V • wV = 0) as required, by equations of MHD. It is now clear that the topological 
properties of the magnetic flux lines can be connected with the topological properties of the velocity fleld lines. 

In this paper we want to understand if it is at all possible to generate a topologically non-trivial distribution of 
hypercharge fluctuations from a topologically trivial collection of the same fields. More physically we are going to 
investigate how to create knots and twists in a stochastic collection of disconnected magnetic flux loops. We would 
also like to understand if this kind of mechanism can or cannot be implemented in some reasonably motivated 
cosmological framework. 

A way of twisting the topology of a collection of independent (hypermagnetic) flux loops, is to have dynamical 
pseudoscalar fields at finite fermionic density. The reason is twofold. On one hand dynamical pseudoscalars fields are 
coupled to gauge fields through the anomaly and, as a result, they induce a coupling between the two polarizations 
of photons and/or hyperphotons. On the other hand, unlike the ordinary electromagnetic fields, the coupling of the 
hypercharge fleld to fermions is chiral ||g] . The generation of hypermagnetic knots can influence the evolution of the 
chemical potential leading to the production of the B AU H . The B AU can be produced provided the electroweak 
phase transition (EWPT) is strongly first order, provided the typical scale of the HK exceeds the diffusivity scale 
associated with the flnite value of the hyperconductivity, and provided the rate of the slowest processes in the plasma 
is sufficiently high. Unfortunately, in the minimal standard model (MSM) the EWPT cannot be strongly flrst order 
for Higgs boson masses larger than the W boson mass ll^]. It is true that the presence of a (topologically trivial) 
hypermagnetic background [Q|9| can modify the phase diagram but cannot make the EWPT strongly flrst order. 
We will argue that our scenario can have some chances in the context of the minimal supersymmetric standard model 
(MSSM) where, incidentally, the rate of the right-electron chirality flip processes can also be larger than in the MSM. 

Dynamical pseudoscalars have been used in the past in order to generate seeds for the magnetohydrodynamical 
evolution |g^,|l|. Our concern here is different and it is mainly connected with the BAU. Suppose that at some 
moment in the life of the Universe dynamical pseudoscalars fields are present. Then, at flnite conductivity and flnite 
fermionic density (i.e. flnite chemical potential) both, the chemical potential and the pseudoscalar particles, can be 
coupled to the anomaly. 



In flat space, the interactions between a pseudoscalar field and an Abelian gauge fleld can be parameterized as 

c^y'^py"' (1-2) 

where c is the pseudoscalar coupling constant, M the typical energy scale and yap = d[ayi3] is the gauge field 
strength expressed in terms of the vector potentials. In our analysis we will be mainly concerned with the case of the 
hypercharge field and with the case of the ordinary electromagnetic field. Imagine now that the pseudoscalar field 
is uniform but time dependent. If our initial state is a (topologically trivial) stochastic distribution of gauge field 
fluctuations the effect of the interaction described in Eq. ( |l.2| ) will be to couple the two (transverse) polarizations 
of the gauge flelds. Since the evolution equations of the two (independent) polarizations are different, there are no 
reasons to expect that the topological structure of the initial gauge field distribution will be left unchanged by the 
presence of pseudoscalar interactions. The inclusion of a (macroscopic) ohmic conductivity will automatically select 
a preferred reference frame (the so called "plasma frame" ) where the flux and helicity associated with the magnetic 
component of the the Abelian background are both (approximately) conserved if the resistivity is small. 

In the plasma frame, the initial stochastic distribution of gauge field fiuctuations can be viewed as a collection of 
unknotted magnetic flux tubes (closed by transversality) which evolve independently without breaking or intersecting 
each others. The inclusion of pseudoscalar interactions will then alter the topological properties of the magnetic flux 
lines by introducing links between independent loops and by also twisting a single loop. Part of the results reported 
in this paper have been recently presented (in a more compact form) in 122]. 

The plan of our paper is then the following. In Section II we introduce the basic description of hypermagnctic fields 
coupled to dynamical pseudoscalars in an empty (but expanding) space. We will then solve the time evolution of the 
field operators in the Heisenberg representation. In Section III we will formally address the problem of emission of 
hypermagnctic knots induced by dynamical pseudoscalars. In Section IV we will concentrate our attention on the 
hypermagnctic helicity production during an inflationary phase whereas in Section V we will discuss the radiation 
dominated phase. In Section VI we will connect the production of (mean) hypermagnctic helicity to the BAU and in 
Section VII we will investigate the regions of parameter space where the generate BAU can be phenomenologically 
relevant. In Section VIII we will examine, for sake of completeness, the problem of global parity breaking induced by 
axionic particles. Section IX contains our concluding remarks. In the Appendix we collected some useful technical 
results. 

II. HYPERMAGNETIC FIELDS AND CHERN-SIMONS WAVES IN EMPTY SPACE 

The Abelian nature of the hypercharge field does not imply that the hypermagnctic fiux lines should have a 
trivial topological structure. We will name topologically trivial the configurations whose field lines are closed (by 
transversality of the field) and unknotted. Conversely, in the gauge 3^o = 0, V • 3? = 0, an example of topologically 
non-trivial configuration of the hypercharge field is the Chern-Simons wave |&t13 



yx{z,t) ^y{t) sin koz, 
yy{z,t) = y{t) cos koz, 
y.iz,t) = 0. (2.1) 

This particular configuration is not homogeneous but it describes a hypermagnctic knot with homogeneous helicity 
and Chern-Simons number density 

nY-^xnY = kQH^{t), 

ncs = -4r^'HY-y^Tr^n\t), 2.2 

where Tiy = V x 3^, T-L{t) = fco3^(i); g' is the U{1)y couphng. Other examples can be found. For instance, it is 
possible to construct hypermagnetic knot configurations with finite energy and helicity which are localized in space 
within a typical size Lg. These examples are reported in Appendix A. 

Suppose now that the Universe can be described, at some epoch below the Planck scale, by a homogeneous, isotropic 
and spatially flat Friedmann-Robertson- Walker (FRW) metric with line element 

ds"^ = g^udx^dx" , gi_,^ = a^{T)T]^^u (2.3) 



(r is the conformal time coordinate and r^^i^ is the ordinary flat metric with signature [+, — , — , — ]). Let us assume 
that dynamical pseudoscalar particles are evolving in the background geometry given by Eq. (^_^). The pseudoscalars 
are not a source of the background (i.e. they do not affect the time evolution of the scale factor) but, nonetheless, 
they evolve according to their specific dynamics and can excite other degrees of freedom. 

The effective action describing the interaction of a dynamical pseudoscalars with hypercharge fields can be written, 
in curved space as 



d^x. 



\g^^d^i^dp^ - 1/(V) - \Yo.pY'^f' + c^Y^pY'^P 



(2.4) 



This action is quite generic. In the case F(V') = (to^/2)?/'^ Eq. (2.4) is nothing but the curved space generalization 



of the model usually employed in direct searches of axionic particles [g3[. The constant in front of the anomaly 
is a model-dependent factor. For e xam ple, in the case of axionic particles , for large temperatures T > mw, the 
Abelian gauge fields present in Eq. (2^) will be hypercharge fields and c = c^ya'/(27r) where a' = g''^ /4:Tr and c^y 
is a numerical factor of order 1 which can be computed (in a specific axion scenario) by knowing the Peccei-Quinn 
charges of all the fermions present in the model pj] . For small temperatures T < raw we have that the Abelian fields 



present in the action (2.4) will coincide with ordinary electromagnetic fields and c = c^-yacm/27r where aom is the fine 



structure constant and c^^ is again a numerical factor. 

The coupled system of equations describing the evolution of the pseudoscalars and of the Abelian gauge fields can 
be easily derived by varying the action with respect to ip and Y"^ , 



D'f/)- 



dV_ 



4M 



Y aY"^ 



V^J^" 



M 



v.^y 



l^y 



v^y'^" 



(2.5) 



where, 



v^y^'' = 



--d„ 



-gY^^^ 



v^y^-^ = 



--du 



^Y 



^ij 



are the usual covariant derivatives defined from the background metric of Eq. 
[-g]'~^'^^a[^f^g°''^^p1j}]. Recall now that 



( |2.3| ) and U^ 



Y 



Oi 



o- £i, YQi ^ a Bi, Yi„ 



-a eijkBk, Yij — a eijkSk 



(2.6) 

(2.7) 



where £i and Bj a rc the flat spa ce fi elds (th e con travariant components obtained by raising the indices with the metric 
given in Eq. (^.3[)). Using Eq. (2.7), Eqs. (|2.5|) can be written in terms of the physical gauge fields 



V'" + 2H^p' - 
V ■ By = 0, 



,dV 



1 



a'' M 

0, V -Ey = 



V y. Ey+B'y 



M 



a 
Vip- By 



V X By ^E'y- — 



iP'By +Vip X Ey 



(2.8) 



(notice that the prime denotes differentiation with respect to the conformal time ). The rescaled fields Ey and By 

are related to £y and Hy as Ey = a^fy, By = o^By- 

We want now to study the amplification of gauge field fluctuations induced by the time evolution of tp. Then, the 
evolution equation for the hypermagnetic fluctuations Hy can be obtained by linearizing Eqs. (2.8). We will assume 



that any background gauge field is absent. In the linearisation procedure we will also assume that the pseudoscalar 
field can be treated as completely homogeneous (i.e.jVV'l ^ V'')- This seems to be natural if, prior to the radiation 
dominated epoch, an inflationary phase diluted the gradients of the pseudoscalar. 

In this approximation, the result of the linearization can be simply written in terms of the vector potentials in the 
gauge yo = and V • y = 0: 



y" - v^y + -rTV^'v X y = 0, 

IVi 

oip a a 



(2.9) 
(2.10) 



where the over-dot denotes differentiation with respect to the cosmic time coordinate t [i.e. a{T)dT = din By 
combining the evolution equations for the gauge fields we can find a decoupled evolution equation for Hy = V x y, 



and we can also decompose the vector potentials in Fourier integrals 



Yix,r)=Y^ 



^W) 



Ya,{k,T)ea,e"'-'' +Y*{k,T)ea,e 



-ik-x 



(2.11) 



(2.12) 



where a = 1, 2 and runs over the two (real) linear polarization vectors e^ whose direction depends on the propagations 
direction k. According to this decompositi on, e i, 62 and 6 3 — k/\k\ form a set of three mutually orthonormal unit 
vectors (i.e. ei x 6*2 = 63). Inserting Eq. (2.12) into Eq. (2.9) we get a system of (still coupled) equations for the 
Fourier components 



r/' + evi - ik^Y2 - 0, n" + ^¥2 + ik^Yi =. 0, 

M M 



(2.13) 



In order to decouple this system of equations we can pass from linear polarization vectors to (complex) circular 
plarizations by defining: 



F±(fc,T)=Yi(fc,r)±zy2(fc,T), Yl{k,T)^Y{{k,T)T^Y2*{k,T) 



In terms of these linear combinations the decomposition of (2.12) becomes 



Y{x,t) 



1 



V2 



E 



d^k 



* ^—ik-x 



where now /? = +,— labels the two circular polarizations 

ei - ie2 



Yp{k,T)epe''-- + Y*{k,T)e0*e-''- 



e+ 



V2 ' 



ei + ie2 



(2.14) 



(2.15) 



(2.16) 



satisfying the usual identities 



e+ • e- = 1, £*+ =£*-, el = e+, 
e+ • eV = 0, e_ • e_ = 0, 

k X e+ ^ i\k\e+, k x e- = —i\k\e^. 



In terms of Y± Eqs. (2.13) can be written as 



YI+ujIY±=Q, Y*1+ojIY1^0, ^l = k^Tk—i'' 



(2.17) 



(2.18) 



Notice that Y± are the normal modes of the action describing the (c lass ical) hypercharge fluctuations excited by the 
evolution of the dynamical pseusoscalar. In other words the action (2^) is not diag onal in terms of Yi^2{x,t) but it 
can be diagona lised in terms of Y±. By inserting the decomposition given in Eq. (2.15) back into the action (p]j), 
by using Eqs. ( p.l7| ) and by integrating over d^x we obtain that the action for the hypercharge fluctuations can be 
expressed as: 



S=fdTL{T), L{t) = f d^kC{k,- 



(2.19) 



^ In order to avoid confusions we denoted with H the Hubble factor is cosmic time and with Ti the Hubble factor in conformal 
time. These quantities cannot be confused with Hy ~ a^Ti-y since, when some ambiguity might arise, we will always Iceep the 
appropriate subscripts. 



with 



£(fc, r) = y[(fc, T)Y*'{k, t) + YL{k, T)Yl'{k, r) - u;lY+{k, r)y;(fc, r) - u;^_Y^{k, T)Y*{k, r) 



(2.20) 



[notice that Y*{k,T) = r+(-fc,r), Y*{k,T) = Y^{-ky)]. 

In order to make more exphcit the meaning of Eqs. ( |2.9| ) and (2.10) let us start from some quahtative considerations. 
Suppose that the potential term appearing in the action (UM is exactly {m? /2)ilP' . Moreover, for sake of simplicity 



we will deal with the case where the co nduc ti vity is completely absent. If i/ ^ ?ti the V' field freely evolves, and, as 
a consequence, we will have, from Eqs. (2.8)-(2.18) 



1 



^'--7 



n'- 



k^^k 



a?M 



1± = 0, 



(2.21) 



and similarly for the other complex conjugate solutions. At H ' 
and typical amplitude ipg: 



ip will start oscillating with typical frequency m 



^{t) 



^0 

0,3/2 



sin [m(t — to)]. 



(2.22) 



Since H < m the effect of the Universe expansion can be neglected during t he os cillating phase as long as we limit 
ourselves to time intervals At<H^^. Moreover, since ip' — atp from Eq. ( 2.18 ) we can see that for ip > u!c/M, 
y+ will be amplified (we defined co ~ k/a as the physical frequency). The maximal amplified frequency will be 
tJinax — V'maxc/Af ~ cm{ijjQ / AI) . In thc same regime Y- is not amplified. Then it is not unreasonable to expect that 
pseudoscalar quantities (like F • V x F, Hy ■ V x Hy etc.) can be produced. Of course, the fact that pseudoscalar 
configurations of the hypercharge fields are produced docs not necessarily mean that they cannot be erased by the 
finite value of the conductivity. We will come back on this point later. 



III. HYPERMAGNETIC HELICITY AMPLIFICATION 



Starting from a stochastic mixture of topologically trivial fluctuations, the time evolution of a generic dynamical 
pseudoscalar excites gauge field configurations which are topologically non-trivial. Suppose that the gauge field 
fiuctuations are in their vacuum state, and suppose that after a period of free-rolling the pseudoscalar field starts 
oscillating at a curvature scale H r^ m. This can happen both in towards the end of an inflationary epoch and/or in 
a radiation dominated regime, as we will show in the next two Sections. 

Particle creation in an external field [E5| plays an important role in the analysis of graviton production due to the 
evolution of the gravitational background [ p6[ and it is also crucial for a rigorous treatment of the normalization of 
density fluctuations in thc context of ordinary inflationary models p^. In quantum optics |2q| similar techniques 
are used in order to describe the photon production in an external (pump) laser fleld interacting with a non-linear 
material. In condensed matter theory these techniques are also employed and their similarities with cosmological 
problems are certainly worth to emphasize [g9| . 

In the Heisenberg description the process of amplification can be understood in terms of unitary transformations (the 
Bogoliubov-Valatin transformations) relating the field operators between two asymptotic vacua. In the Schroedinger 
picture the amplification process is described through the action of a time evolution operator connecting the initial 
vacuum state to the final (multiparticle) "squeezed" vacuum state [g6 28 3^ . 

Th e cla ssical (canonical) Hamiltonian for the gauge field fluctuations can be obtained from thc canonical action of 
Eq. (|^ 



H{t) 



(fkn{k,T), n{k,T) 



7r_7r_ 



'<y;Y+ 



■uj'^Yl 



y_ 



(3.1) 



where 7r± and tt^. are the canonical momenta and where ijj\ are not necessarily positive definite. From Eq. (3.1) the 
evolution equations ( 2.18| ) can be directly obtained in the Hamiltonian formalism. Our problem is analogous to the 
quantization of two decoupled (complex) harmonic oscillators and it is closely analogous to the quantization of the 
two circular polarization of the electromagnetic field in the radiation gauge. The only (but crucial) difference is that 
the two polarizations obey (in our case) two different evolution equations (i.e. U!\ 7^ i^-)- In the limit tp' = (i.e. 
CJ+ = uj^ ) the ordinary "electromagnetic" case is completely recovered |g^ . 
Promoting the classical normal modes to quantum mechanical operators 



T^+{k,T) -^ 7r+(fc,r), Y+{k,T) -^ Y+{k,T), 

^_(fc,T)^^_(fc,T), y_(fc,r)^y_(fc,r), (3.2) 

we can impose the conimutation relations 

[Ypik,T),np,{p,T)]^iSpi,,S'^^\k-p), [Yi3{k,T),Yp,ip,T)]=0, [^^(fc, r), V(p,t)] = (3.3) 

where f3 — +, —. The dynamical evolution follows the Hcisenberg equations of motion given by 

Y;,=z[H,Yp], rr'p=t[H,7Tp], (3.4) 

where H is the Hamiltonian operator. A consequence of these equations is that the operator Y± obey , in the Heisenberg 
representation, the same evolution equation of the related classical quantity derived in Eqs. ( 2.18| ). 
We define the creation and annihilation operators as 



r+(fc,r) 



1 



V"^ 



a+{-k,T) + a\{k,T) 



Y-ik,T) 



1 



-.fhJZ 



7r+(fc,T) = -«y^ 



a+(k,T) -a\_{-k,T) 



7r_(fc,T) =-i\/ — 



a_ (— fc, r) + a_ {k, t) 
a^{k, t) — alL(— fc, t) 



(3.5) 



Notice that a±[—k,'q) are annihilation operators for an excitation with momentum —k whereas a±{k,r]) are annihila- 
tion operators for an excitation with momentum +k. The simultaneous presence of both sets of operators guarantees 
the explicit three momentum conservation of the whole formalism. The commutation relations obeyed by these 
operators are the standard ones, namely 



[a,s{k,T),dl,,{p,T)] ^5pf3'S^^\k-p), [a+(fc,T),a_(p,T)] = 0, 
and the Hamiltonian operator can then be written as 



Hir) = E / 



d^fc 



^l^i'^) [ o^pQ^i T)0'l3{ki t) + a/3(-fc, T)ap{-k, t) + 1 



13=-,+ ' 
The field operators before the amplification took place can then be written as: 



Y, 



0-+,inJin + '^+ in/ii 



, Y-. 



0--,in^in + '^^ in^in 



The same decomposition can be done for the operators after the amplification processes has taken place: 



Y. 



-{-,out 



0'-{',out9out ~r (^J^^outdout 



r_ 



^—.OUt^OUt \ Q._ nii.t^n 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



We can formally write the out-going mode functions in terms of the in-going mode functions : 

goutir) = c+Ui{t) + C^fl^{T), 
Goutiv) = c+F„,{t)+c^F*Jt). 

Since y^'™ and Y^'°^* represent the same solution they must be equal: 
Inserting the explicit form of the solutions we have that 

O'+.inJin ~r ^--i-^inJin ~ '^+,outgout ~r 0,+ .out9out 



(3.10) 
(3.11) 

(3.12) 



Using Eqs. (3.10) in Eq. ( 3.12| ) we obtain that the creation and the destruction operators after the amplification 
(^±,in) are related to the ones before the amplification (a±_i„) by the following Bogoliubov-Valatin transformations: 






(3.13) 
(3.14) 



Once the evolution of the pseudoscalar field "ip is completely specified, the coefficients c± and c± can be explicitly 
computed. It is useful, at this point to re-write the full expression of the hypercharge after the amplification took 
place 



1 



Yout {x, t) ^ —= 



Sk 



a,+,out9out + a^ ^g^^^ K+ + a-,outGout + «- out^mit ^ 



V2J 7(2^ 

^t . ,~ i^ ,/^t r* .- r 1^ 

^-\-,out9out ~^ ^-\-, out 9 out \^- ~^ \ ^-.out^out ~^ ^-,out^out K+ 



Ak-x 



^—ik-: 



(3.15) 



We have now all the ingredients in order to estimate the expectation values of the pseudoscalar quantities we are 
interested in. As it was previously stressed EflJllt] a good gauge-invariant measure of global parity violation is 
represented by Hy ■ V x Hy- This operator emerges naturally in the study of the evolution of hypermagnetic fields at 
finite conductivity and finite fermionic density. As we will show in the next two sections, the same provides a good 
gauge- invariant measure of parity breaking both for T >Tc a nd fo r T < Tc- 

Then we can compute our expectation value by using Eq. (3.15) and the Bogoliubov transformations (in order to 
relate incoming and outgoing modes): 



\Hy -V X Hy\0) 



d^P 3 

(27r)3^ 



(|c+|^ + \c-\')Mr)f:jT) - (|5+|^ + |5_nF.„(r)i^;(r) 



(3.16) 



where |0) denotes the vacuum state of the hypercharge field. Suppose that our initial field configuration is given by 
a stochastic background or by a thermal mixture [p3 . Then /i„ (r) and Fin (t) will be 



fUr) = ^^~''^ 



FUr) = :^e-(^-+^). 



(3.17) 



If there is no initial phase shift (as it should be) (/S = and then the two mode functions are equal. Thus, going 
to Eq. ( 3.16 ) we have that the whole expression is zero provided c± ^ c±. Indeed we showed that because y+ and 
y_ obey different classical equation the Bogoliubov coefficients with and without tilde will be in general different. 
The specific form of the Bogoliubov coefhcients will of course depend upon the specific cosmological model, upon the 
specific particle physics candidate for -0. However, even in the general case we can carry on the calculation one step 
further. The amplification coefficients are represented, in our formalism, by c_ and c_. Thus, the logarithmic energy 
spectrum of the hypermagnetic helicity can be written as 



dJiey 
dloguj 



1 



-OJ^g{u:,T)\A{uo)\\ 



with 



Qi^.r) 



|c_(a-)p - |5_(a-)p - c^c^flir) + c+c^ft^ir) - ^.c^Flir) - £+c1F*,'{t) 



(3.18) 



(3.19) 



where we defined, as usual, the physical frequency lo = k/a. Notice that in deriving Eq. (B.18) we used the relation 



|c+| 



1 (and analogously for the tilded quantities) which is a trivial consequence of the unitarity. 



By defining the logarithmic energy spectrum of the initial distribution of hypemagnetic fields ||11| 



p{uj) 



dp\ 



dlogw 






Eq. (3.18) becomes 



dJiey 
dlogw 



= ujQ{u!,t) p{lu). 



(3.20) 



(3.21) 



Even if the amplification related to the dynamical pseudoscalars is of order of 1 , this last expression tells us th at the 



obtained helicity can be large provided the initial state is not the vacuum. We will see that the result of Eq. ( 3.21 ) 
can be quite relevant for our applications. 

Therefore, in this section we derived a number of results. As previously noted |20[| we obtained that dynamical 
pseudoscalars induce an interaction between the polarizations of Abelian gauge fields. This coupling has the effect of 
modifying the evolution of the two circular polarizations in such a way that the magnetic helicity (initially zero) can 
be dynamically generated. 

Up to now our results do not take into account the possible effect, in the final state, of some effective Ohmic current 
which can account (macroscopically) of the dissipative effects associated with a plasma at finite conductivity. This 
will be an important point in the next two sections. On one hand the finite value of the conductivity provides a useful 
(and physicalUy motivated) ultraviolet cut-off for our helicity and energy spectra. On the other hand we will stress 
that, also at finite conductivity, the Chern-Simons density can be directly related to the helicity. 

IV. HYPERMAGNETIC KNOTS PRODUCTION DURING AN INFLATIONARY PHASE 

The ways our mechanism can be implemented depend upon the early history of the Universe. We will firstly consider 
a simplified model where an inflationary phase occurs H. The inflationary phase terminates at the time ri and the 
ordinary (radiation-dominated) phase settles in. 

Notice that ip can oscillate both in the inflationary phase and in the radiation dominated phase. These two logical 
possibilities involve two different physical pictures. In the case of inflationary production of hypermagnetic helicity 
the oscillations will only be damped by the Universe expansion. If, on the contrary, ip oscillates during the radiation 
dominated phase we can expect that the evolution of the hypercharge fields will be damped by the finite value of the 
conductivity. In the present and in the following Sections we will analyse these two regimes. 

Suppose that our (massive) pseudoscalar field evolves during an inflationary phase following its equation of motion 
which we write (in cosmic) time as 

V^ + Si/V- + mV = 0, (4.1) 

where H is (approximately) constant and it corresponds to the maximal (curvature) scale reached during inflation. In 
this model we have essentially three parameters. One is Hi/Mp fixing the maximal curvature scale during inflation. 
From the contribution of gravitational waves to the Cosmic Microwave Background (CMB) anisotropy we know that 

ll^^lO-, (4.2) 

where Hi is the curvature scale at the end of the inflationary phase. Moreover, we also know that the EWPT takes 
place when T > 100 GeV. Since, by that time, the Universe was dominated by radiation we also have to assume that 

^ > 10-''. (4.3) 

Mp ^ ' 

The second parameter of our model is m/ Hi which essentially fixes the amount of amplification of the hypercharge. In 



Eq. (2.18) the instability develops for k '^ fcmax ~ c{'ipo/M)ma. The amplification occurring for k ~ fcmax in the mode 
function will be of the order of Y^{k,T) ~ exp [c{'ipo/M)maAT] where aAr ^ At ^ H^^ is the typical amplification 
time. The typical growth of y+ will then be proportional to exp [c{ipo/M){m/Hi)]. As we can see the amplification 
depends also upon (-tpo/AI) which is the third parameter of our analysis. The coupling constant c appearing in the 



action (2.4), is in principle, determined by the specific particle physics model. We want immediately to notice that the 
amplification occurs for m^ Hi. At the same time the coherence length of the amplification processes is determined 
by iinax ^ ^mix ~ c""^ (M/^/'o)m~-^ and, therefore, the larger the amplification is the shorter is the coherence length 
of the amplified field in horizon units, lim > Hi effect of the Universe expansion can be neglected at least in the first 
approximation. 

This qualitative analysis can be refined by a numerical study which includes the Universe expansion. For this 
purpose we can write the hypercharge evolution in cosmic time. The system which should be solved numerically is, 
in general, 

c • H^ H' 



y± 



M^ 4 2 



2/± = 0, 



*4 

y± 



:H' - Ih 



4 2 



* = 0, 



(4.4) 



In our specific case this second order differential system can be simplified since H = 0. We studied numerically this 
set of equations by reducing it to a first order linear differential system. Some examples of our results are reported in 
Fig. y. The crucial point in our analysis is the dete rmination of the maximally ampplified mode Wmax- Once this is 
obtained we can solve the system given in Eq. (4.4) to include the Universe expansion. In the approximation where 
m^ H {i. e. expansion negligible compared to the oscillations) we have that the Eqs. (4.4) can be further simplified 
since ip ~ a~^^^ sin [m{t — to)]. We get a Mathieu-type equation ^5 36 



fy± 



S =F ecos2z 



y± = o, 



(4.5) 



where 5 — 
k = 0,1,2., 



(2cj/m)2, e = citpo/M)(2Lo/m) and z = (mi)/2. Note that when e = and S = {2k + if (where 
) there is solution of Eq. ( |4.5| ) with period 27r but when e = and 6 — (2fc)^ there is a solution of Eq. 
(4.5) with period it. Therefore, in the plane {5, e) the instability bou ndaries will cross the e = at the points defined 
by ^ = {2ky or S = (2fc + 1)^. By studying the instability of Eq. (15) in the case of e ^ (but still e < 1 ) we get that 
the first instability region (corresponding to (5 ~ 1) occurs for a region in the plane (J, e) bounded by (5 — 1 ~ ±e/2. 
Inserting the values of S and e we get that, for small e this region is centered around Wmax ^ {c/2){tpo/M)m. Of 
course there will be also other bands defined by their intersections on the e = axis {5 = 4, 9, 16, 20...). Those bands 
correspond to higher frequencies but they are narrower than the first one defined by the 6—1 intersection [BSpq] . 
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FIG. 1. We report two numerical solutions of the system given in Eqs. (4.4). We assume uj = cjmax- At the left we report 
the case of c = 0.01 and H/m = 0.001. At the right we report the case with c = 0.001 and H/m ~ 0.00013. For the two cases 
we would get, from our considerations following Eq. ([4.5|), amplifications of the order of e^ (left) and e* (right). The damping 
due to expansion settles in quickly (roughly after H~ in mass units). The agreement our our estimates with the numerical 
evaluation gets worse for frequencies far from the maximal amplified frequency. 

According to this analysis we can say that Wmax leads to an amplification of the order of exp [{c/2){iljo/M)mAt], 
where At ^ H~'^ . As noted in previous studies of a similar equation, the Universe expansion will soon shut-off the 
instability occurring in the e < 1 case |34 . 

In Fig. 1^ we report the solution of the system given by Eq. (14) where we took lj — uj-max- The Universe expansion 
has been included. If we take the case with c = 0.01 and m/H = 1000 we get that the amplification should be 
e^ — 148. Numerically we get that during one Hubble time the maximal amplification is roughly 150 (left picture in 
Fig. M. However the amplification is quickly shadowed by the Universe expansion which becomes effective roughly 
after a time of the order of H~^ in mass units. 

Thus pseudoscalar oscillations are quite effective in changing the topology of a stochastic background by inducing 
correlations between the two independent polarizations of the gauge field. However, they are not effective in inflating 
the amplitude of the background. Therefore, as we will see even better in the case of oscillations occurring in a 
relativistic plasma, one of the crucial assumption of our considerations will be the existence of a primordial stochastic 
background of hypcrmagnetic fields whose topology can be eventually twisted by the pseudoscalar oscillations. 

In Fig. n^ we completely ignored the role played by the finite value of the conductivity. During the infiationary 
dominated regime the effective conductivity is equal to zero since no charged particles are present in the Universe. 
When the Universe reheats charged particles are generated and, therefore, an Ohmic current is developed. This effect 
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can be modeled by thinking about the situation where the conductivity is zero prior to the onset of the radiation 
dominated phase but it gets larger and larger as soon as charged particles are created. 

Needless to say that the the smooth transition of the conductivity from zero to a finite (large) value is quite 
important for our considerations. Suppose, as a warm-up, that we are in flat space and suppose that the conductivity 
jumps from zero to a finite value at some time i*. Let us also assume, for simplicity, that there are no pseudoscalar 
interactions in the game. Then the evolution equation of the hypercharge fields will be the usual D'Alcmbert equation 
for t < ti, but it will be modified as: 



y" + <ycy' - v^:^ = o. 



(4.6) 



for t > t^,. If we match the two solution at i = t^, we simpl y g et that the hypermagnetic fields are left unchanged by 
the transition whereas the hyperelectric fields are damped |43] : 



n>{u) = n<{h), £>iu) =. e-"^**£<(t.). 



(4.7) 



where the subscript > (<) denotes the value of the fields after (before) the transition occurring at i*. Since pseudoscalar 
interactions are absent in this example we also have that (3^+ — 3^_) will always be zero. 

From a physical point of view we can expect that the effect of the finite value of the conductivity will be to fix a 
definite sign for (>+ — Y^). By ignoring the effect of the conductivity we can argue (see Fig. nl) that (y+ — Y-) is 
certainly amplified but, at the same time, the sign of this quantity is not well defined. In order to clarify this point 
let us study the full system of equations at finite conductivity 



n' 



aY^ 



k^Tk-^' 



Y±=0, 



ip" + 2nij' + m^a^^ = 0, 



(4.8) 



where we remind that the conductivity, in curved space, evolves as (j(t) = aca{T). In order to model a smooth 
transition from the inflationary epoch to the radiation dominated epoch let us assume that the geometry evolves 
continuously from a de Sitter (or quasi do Sitter) stage to a radiation dominated stage of expansion. This can be 
achieved with a specific choice of the scale factor, namely 

(4.9) 

Notice that for r -^ +oo we have that a{T) ~ r. For r -^ — oo we have that a{T) ^ t^^. For r — > we have that 
a{T) '^ constant. So we can say that this geometry evolves smoothly from a de Sitter phase to a radiation dominated 
phase. The Hubble parameter in conformal time can be easily computed and it is also a smooth function of r 




nir) 



1 



yA^T' 



(4.10) 



The duration of the transition period between the de Sitter phase and the radiation dominated phase is controled 
by Ti. By increasing ri the transition period gets longer. By decreasing ri the transition period gets shorter. This 
aspect can be appreciated by looking at the time evolution of H(r) (see Fig. g). 
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FIG. 2. We illustrate the Hubble parameter TI = a' /a computed from Eq. (4.9). By increasing ri the transition regime gets 
longer. For ri = 1 the Hubble parameter is reported with the full (black) line. For larger values of ri (i.e. ri = 5 [ thick dashed 
line] and ti = 10 [full thin line]) the height of the Hubble parameter decreases and the transition regime gets broader. 
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The behavior of the conductivity follows from the evolution of the geometry. Therefore, in this model, we will have 
that the conductivity is vanishingly small during the de Sitter phase and it grows linearly in conformal time during 
the radiation dominated phase. By rescaling with respect to the mass in order to have fully dimensionless quantities 
we can easily re- write Eqs. (4.S) in the form of a first order differential system 



dq 



dq- , ■, d.u 



dp+ f(Tc\ I ^ 

dx V "^ / 


^'-^i^-i-)_ 


9+ (a;), 


— — = - — a(x)p_ - 
dx \m J 




9- (a;). 


dv 2 v{x) 2t \ f ^ 
dx ^x^ + x\ 


Vj_ and u — -ip and v — ■(/;'. 


Notice, moreover th£ 



(4.11) 

where we denoted q± = y±, p± = Y^ and u — -ip and v — ^'. Notice, moreover that x = mr, x\ = mT\ and e = kjm. 
We can observe that there are three important parameters for our discussion. The first one is, Tii = 'H(t\) ^ Tj~^ 
the Hubble parameter at the end of the inflationary phase. The Hubble parameter in cosmic time is simply related 
to Til since Tii = a\H\ '^ (1 + \/2)Hi. Moreover, as we noticed before, in order to be compatible with the large 
scale observations in the microwave sky we have to require Hi/Mp < 10^*^. Notice that Hi always appears in the 
combination xi = mri ~ m/Hi. Consequently, if we w ant la rge pseudoscalar oscillations prior to ri we have to 
require Xi <^ 1. The last parameter appearing in Eqs. (4.11) is ac/m. Numerically CTc — fo^o where CTq — 70- 
100 and Tq is the temperature of the Universe at the beginning of the radiation dominated epoch. Notice that 
ac/m = {ac/Hi){Hi/m). Therefore, taking into account that Hi < 10~^Mp and that we want Tg as high as 10^° 
GeV we have to require that ac/m < 1. Suppose, for instance, that Hi ~ lO^^Afp, Tq ^ 10^° GeV, m/Hi = 100. 
Then, we get that ac/m = cro/100 ^ 0.7. We will keep xi and ac/m as free parameters in the integration. Notice 
finally that we will focus our attention on the maximally amplified modes and, therefore, we will fix e ^ c/2. 
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FIG. 3. We illustrate the numerical integration of the system given in Eqs. ( 1.11 ). At the left we report the cases where 
ac/m = 0.1 and xi = 100 (full thick line) and xi — 10* (dot-dashed Une). At the right we fix xi — 100 and we choose, 
respectively, ac/m = 0.1 (full thick line), ac/m = 0.01 (full thin line) and ac/m — 10~* (dot-dashed line). 

In Fig. and H we report the results of the numerical integration for different sets of parameters. In Fig. I| (left 
plot) we fix ac/m, ~ 0.1 and we choose xi = 100 (full line), xi = 10"*. We see that (F+ — F_) acquires a definite 
sign as we expected. Moreover the transition regime gets broader as soon as we increase xi. This can be understood 
by recalling that by increasing ti the transition regime gets more significant. What happens if we change ac/m? In 
Fig. g (left plot) we fix xi = 100 and we decrease the value of ac/m. The the full (thick) line corresponds to the 
case ac/m = 0.1 , whereas the dot-dashed line corresponds to the case ac/m = 10^*. The full (thin) line corresponds 
to the case ac/m. = 0.01. We clearly see that in the limit ac/m, -^ we recover partially the case of Fig. |l| since 
(yi)_ — Y-) does not have definite sign when the radiation dominated phase settles in. The fact that ac/m cannot 
be too small can be understood physically by recalling that the small ac/m limit is equivalent to a small reheating 
temperature. What happens when ac/m, > 1. We can expect that if the conductivity is too large at the beginning 
of the radiation dominated epoch the damping term will be dominant and the amplification will be reduced. This 
aspect is illustrated in Fig. S. 
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FIG. 4. We illustrate the integration of the system given in Eqs. (4.11) in the case where xi — 100, ac/m — 1 (full thin line) 
and in the case where Xi = 100, (Jc/m = 5 (dot-dashed line). 



The limit (Tc/'tI' S> 1 simply means that we are in the damped regime which will be specifically addressed in Sections 
V and VII. 

In the sudden approximation, the Bogoliubov coefficients describing the amplification of the mode function when 
the background changes from the inflationary phase to the radiation dominate phase can be analytically computed 
in the framework of our model. When the modes of the fields evolve in the radiation dominated phase we will also 
have to take into account the effect of the conductivity and this will be the second step of our calculation. We will 
solve the evolution of the mode function for modes close to Wmax. For t < tq the pseudoscalar does not oscillate but 
when its mass is of the order of Hi the inflationary oscillations will begin until the radiation dominated phase starts 
at r ~ Ti . The mode function corresponding to y+ for k > fcmax evolves in the three regions as 



1 



'''^^^ = Tk 
'"'^^^ ^ Tk 



, r < To 



c+e-'''^ + c_,e'''^ 



To < T < Tl, 
T > Tl, 



(4.12) 



where /3 = c/2(il}o/M)ma. We want to stress that the solution for To < t < ti holds for modes k ^ /cmax and for times 
aAT<H~^ (where At ^ ti — To) but it should not be trusted outside of this range as discussed in Fig. ^. Notice that 
c± are exactly the Bogoliubov coefficients we are looking for since they relate the incoming mode function with the 
outgoing mode function. They can be determined by matching the three solutions (and their first derivatives) in the 
three regions defined by tq and ti . This can be easily done and the result is 



1 



■ikAr 



2 cosh [/3At] 



^-f)sinh[^AT] 



'- = -r 



-ikAr 



l + f )sinh[/3AT] 



(4.13) 



where At = t^ — to defines the duration of the oscillatory phase. During the time the 1+ is amplified we also have 
that the Y- will oscillate and the corresponding Bogoliubov coefficients will not describe amplification but mixing 
between the positive and negative frequencies. They are 



1 



ikAr 



2 COS [/5At] + i 



'ikAr 



(3 



sin [/5At] 



sin [(3 At] 



(4.14) 



Our expressions of the Bogoliubov coefficients satisfy |c+ 1 



1, and jc-)- 



;; |2_ 



15-P 



1 as required by the unitarity 



of the Bogoliubov transformation which preserves the commutation relations between the creation and annihilation 
operators in the asy mptot ic regions. 

According to Eq. ( [3.1^ ), the mean helicity spectrum in the case A{k) = 1 will be: 



dHey 
dlogu; 



'^max c(m/Hi)(i)o/M) 
2 



(4.15) 
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Thus according to Eq. ( 3.2l| ) we will have: 



r('^max)''|-4(Wmax)|' 



(4.16) 



In Eq. (4.15) we assumed that the hypermagnetic field modes are re-entering in the radiation dominated epoch but 
we did not take into account properly the damping effect of the conductivity. If the field modes are amplified (as in 
our case) during a phase where the conductivity is strictly zero (since no charged particles are present) the only effect 
of the conductivity (arising when the field modes enter in the radiation dominated phase) is to cut-off (exponentially) 
all the momenta larger than the typical momentum associated with the conductivity. However, if the pseudoscalar 
oscillations occur in the radiation dominated phase the effect of the conductivity can be even more important, as we 
will show in the next section. 



V. HYPERMAGNETIC KNOTS PRODUCTION DURING A RADIATION PHASE 



If pseudoscalar particles are present prior to the onset of the electroweak epoch they can offer a mechanism for the 
production of Chern-Simons condensates whose decay can seed the BAU. An example of such a particle is the axial 
Higgs whose mass enters, after supersymmetry breaking, in the expression of the neutral Higgs boson mass ||3l| , |32| . 
Inspired by the axial Higgs we can require the pseudoscalar mass to be larger than 300 GeV q If m > 300 GeV the 
MSSM Higgs sector is not too different from the one of the MSM. 

We want to extend the discussion of hypermagnetic helicity generation to the case of a radiation dominated epoch. 
Suppose now that at some higher curvature scale a stochastic distribution of hypermagnetic field has been created. 
Then, suppose that at some lower energy scale the pseudoscalar oscillations will start developing. In this case the 
evolution of the hypercharge fields is 



Y^ 



aYi 



k =F kc 



Y± = 0. 



By eliminating the first derivative through a field rescaling we get 



y'i 



A-r 4 2 



1 r 

y± = 0, Y±= y± exp [-- / a{T')dT']. 



(5.1) 



(5.2) 



If an initial (stochastic) hypermagnetic field is present, then the pseudoscalar interactions will be able to twist the 
topologically trivial distribution. If no initial stochastic distribution is present, the obtained hypermagnetic helicity 
will be quite small. In other words we find that at finite conductivity the pseudoscalar interactions are only efficient in 
twisting a typologically trivial hypermagnetic distribution. Indeed, as we previously discussed, we can always define 
a conductivity scale at any time in the radiation dominated epoch. All the modes of the fields with typical frequency 
larger than the typic al co nductivity frequency are essentially washed away whereas the other have some chance of 
being amplified. Eq. (5_^) becomes, in the large conductivity limit. 



aYi 



k =F kc 



M 



Y±^0, 



and the corresponding solution can be written as 



Y± 



^dr' 



^TI ^-M4''ir')dT' 



A{u;) 



(5.3) 



(5.4) 



where the initial amplitude A{lo) coincides with the typical amplitude of the stochastic hy perm agnetic distribution. 
The maximal value of the hypermagnetic helicity in critical units can be obtained from Eq. (5.4) with the result that. 



(iJy ■ V X Hy) 
crpc 



2ao M mJ^^'^'^'''' 



(5.5) 



^Assuming no invisible decays, the direct limit is m > 65 GeV (at 95% confidence level and in the case tan/3 > 1) |33 
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where ctq = (Jc/Tc and a = UcCi and A^ is the increment of ip from its initial value. Typically (without fine-tunin 



g 



the initial amplitude of tp at the beginning of the oscillations) we will have that A-i/; ~ M . Notice that in Eq. (5.5) 
we also defined 

Pc 1^ Pc ^a -'''-' -'Wo 

as the critical fraction of energy density present in the form of a stochastic background of hypermagnetic fields for 
uj ^ tJmax- Eq- (pj) holds for frequencies smaller than typical conductivity frequency 



^oir) ~ ^a^^/Nlff^l^T (5.7) 



M 



p 



It is quite relevant to notice that in Eq. (5.5) the exponential damping holds for all the modes uj > lo„. The only 
way of getting a huge amplification (at finite conductivity) would be to require Wmax ^ ^a- However, in this limit, 
the conductivity erases all the field modes. Therefore we are forced to discuss the case Wmax < i^o-- This limitation on 
the amplified physical scales applies for all the temperatures (in the symmetric phase of the electroweak theory) where 
helicity generation takes place. 

If the conductivity is finite it does not make sense to consider the case of amplification of the hypermagnetic helicity 
from tha vacuum fluctuations. The reason is that, in this case, ^'(winax) is extremely minute. This can be easily seen 
in general. Suppose, in fact that Wmax can be made as large as O.lwo-- Then we get that the vacuum fluctuations will 
be in this case 

10-%/ T y 

r(u;.„ax)'^^^(^^j (5.8) 

which is very small for in the symmetric phase of the electroweak theory and in the TeV range. 

VI. BAU FROM HYPERMAGNETIC HELICITY 

We are interested in the possible effects of the hypermagnetic fields on the electroweak scale and then we will 
mainly consider the plasma effects associated with the dynamics of Abelian gauge fields in the symmetric phase of 
the electroweak theory. This is consistent with the assumptions made in the previous section where we focused our 
attention on hypermagnetic fields amplified and/or re-entered during a radiation dominated phase. 

The description of the plasma effects arising in the symmetric phase of the electroweak theory has to rely on some 
generalization of the MHD equations. In pd]] it was argued that such a generalization can be provided by the equations 
of anomalous magnetohydrodynamics (AMHD) where the effect of finite conductivity and finite chemical potential 
can be simultaneously described. In the present analysis we have to generalize AMHD to the case of pseudoscalar 
interactions and non vanishing bulk velocity of the plasma. Our variables are the hypermagnetic and hyperelectric 
fields, the right electron chemical potential and the pseudoscalar field ip. Of course these equations have to be 
supplemented by the evolution equation for the bulk velocity of the plasma. The complete set of equations reads 

[{p + p)v\' + v-V[{p + p)v] +vV ■[{p + p)v] = -Vp + Jy X Hy + yylV^w + -V(V ■ v)], ri = {p + p)v 

O 

E'y + Jy = ^PRaHy + —WHy + V^ X £;y] + V X Hy, (6.1) 



H{. = -V X Ey, V-Hy^O, 



M 



2 2; V72 ;. C Ey ■ Hy 



V'^ + SHV^' + m^a^V- VV = -— ^ , cr = aac, 

M a^ 

^+^-[{p + Pn=EyJy, (6.4) 
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where v is the bulk velocity of the plasma, p — a^pr and p = a^pr are the (rescaled) energy and pressure densities of 
the radiation dominated fluid. Notice that v and 1/cr are the thermal and magnetic diffusivity coefficients whereas Dn 
accounts for the dilution of right electrons in the case where the chemical potential is not completely homogeneous. 
Using the fact that Jy '--' V x Hy the Navier-Stokes equation, for scales larger than the thermal diffusivity scale, 
becomes 

[{p + p)v\' + v-V[{p + p)v\+vV- [{p + p)v\ = -V[p + ^\Hy\^] + [Hy ■ V]Hy, (6.5) 

where we used the fact that V x Hy x Hy = ^^^\Hy\^ + [Hy ■ ^]Hy- Notice that we will be always concerned with 
the case where p — p/3 ^ l^yp (i-G- vanishing hypermagnetic pressure). Similarly, the continuity equation will be 
modified as 

^+^-[ip + p)v] = -[^xHY]'. (6.6) 

or a 

Since we are dealing with modes of the hypercharge field with momentum k smaller than 

ka-\I^T, A'fo = ^^ ^' ^ 7.1 X lO^^GeV, (6.7) 

the source term appearing in the continuity equation can be ignored and we can indeed assume that our plasma is 
incompressible (i.e. V • v = 0). In this case, consistency with the continuity equation requires that p and p depend 
only upon t. Then, Eq. ( p.q ) can be further simplified 

dv ^ [Hy ■ V]Hy ,.., 

—- + [v ■y\v ^ — — . (6.8) 

dT ' ' [p + p] 

With the same set of assumptions we can also obtain a generalized version of the hypermagnetic diffusivity equation, 
namely 

— ^ = V X [prHy] - ttV X ^'i/y + V X {vxH) + -V^i/y. 6.9 

where, as in the previous Sections, we focus our attention on the case where IVV'I <C "0'. Again, li k < k^ and if the 
pseudoscalar oscillations already took place, Eq. (|6.9|) can be written as 



BH 

—^ + [v-V]Hy^[Hy-^]v, (6.10) 

OT 

where we re-expressed V x (v x Hy) acc ording to the usual vector identities. Under these approximations a fully 



nonlinear solution of Eqs. ( 3.10 ) and (6^), is given by the Alfven velocity ||l 



Hy Hy , , 

v^± , ^±^ 6.11 

Vp + p Vp 

This result seems to imply that if the hypermagnetic field distribution is tangled (for example as a result of the 
pseudoscalar oscillations), then also the velocity field might be tangled in the same way. The arguments we presented 
up to now hold in the case where Jy x Hy 7^ 0. We want to point out that this is not exactly our case. Indeed, 
as we showed in the previous Sections, the maximally amplified hypercharge modes are k ~ fcmax- For these modes 
{Hy ■ V X Hy) is maximal and, therefore we can also argue that (V x Hy x Hy) ^ 0. Thus, we can argue that 
our system is (approximately) force free |£^. Now, if the hypermagnetic field is not completely homogeneous we 
could argue from the general expression of the Navier-Stokes equation, that the hypermagnetic inhomogeneities could 



drive inhomogeneities in the energy density according to a well known observation by Wassermann |37|. Even if the 
hypermagnetic field would not he force-free we could argue that since the electroweak epoch occurs, in our scenario, 
deep in the radiation epoch, the force term can only seed a decaying mode for the density fluctuations Q. In the 
case of a force free field the Navier-Stokes and diffusivity equations acquire a symmetric form |l5| 
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^^ =Vx{vxHy) + -V^Hy 
OT a 

^ = V X (t/ X cJ) + i/V^w. (6.12) 

OT 

Let us now move to the analysis of the kinetic equation for the chemical potential. We could introduce, in principle, 
a chemical potential for all the perturbative strong and weak processes, Yukawa interactions of leptons and quarks 
which are in thermal equilibrium for temperatures T > T^. However, the processes really crucial for our purposes 
are the slowest perturbative processes related to the C/(l)y anomaly, namely the processes flipping the chirality of 
the right electron which are in thermal equilibrium until sufficiently late because of the smallness of their Yukawa 
coupling 1^ . Therefore the right electron number density will be diluted thanks to perturbative processes (scattering 
of right electrons with the Higgs and gauge bosons and with the top quarks because of their large Yukawa coupling) 
but it will also have a source term coming from the anomaly contribution. Notice that V is the (perturbative) chirality 
changing rate. Moreover, there is a second rate appearing in Eq. ( p.3| ) which is proportional to the hypermagnetic 
energy density. The main effect of Th is to dilute the fermion number even in the absence of perturbative processes 
flipping the chirality of right electrons (i.e. F = 0). 

In Eq. (^.3|) we see that the right electron chemical potential can be diluted also by diffusion. The diffusion constant 



Dji appearing in Eq. (6.3) defines the relation between the diffusion current Jji and the gradient of the right electron 
number density 

Jr = -Dr^ur. (6.13) 

The right electron number dilution equation can then be written, in flat space, as 

duR 5'^ 



dt 47r2 



£Y-nY~T{nR-n'j^)^V-jR. (6.14) 



Using Eqs. (6.1)-(6.2) together with Eq. (6.13) into Eq. (6.14) we can obtain Eq. (6.3). The right electron diffusion 




coefficient can be roughly estimated by taking into account the interactions of right electrons with hyercharge fields 
and it turns out that the typical momentum for which diffusion becomes effective is of the order of 

(6.15) 

For field modes k < ko the diffusion term can be consistently neglected in the study of our problem. The hypermag- 
netic helicity appears to be the source term in the right electrons dilution equation. This is a quite relevant point for 
our considerations since a sizable mean magnetic helicity (possibly generated through psudoscalar oscillations) can 
influence the BAU present in the symmetric phase of the electroweak theory. In our considerations we will assume 
that the hypermagnetic diffusivity and the thermal diffusivity are of the same order. More precisely one can define 
the Prandtl number which usually measures, in the context of ordinary MHD, the relative balance between thermal 
and magnetic diffusivity 

Pr,n = ua='^. (6.16) 

a' 

The approximation of unitary Prandtl number is often employed in 3D MHD simulations [ [l5[ . The approximation of 
unitary Prandtl number is often employed in 3D MHD simulations |15[| . In this limit, hypercharge modes k < k^ are 
not dar npe d by thermal diffusion. In our case (as in the case of realistic MHD simulations) Pr,„ > 1. 

Eq. (|6.3[ ) can be discussed in two physically different regimes. The first one is the case where F > Th, whereas the 
second case is the one where F < Th- In the first case the rate of dilution of the chemical potential is dictated by the 
perturbative processes, whereas, in the second case, it is determined mainly by the Abelian anomaly. 

The preliminary analysis of BAU generation through Chern-Simons waves pT| showed that the case where F > Th 
may lead, potentially, to a higher BAU and therefore we will focus our attention on the case where F > Th commenting, 
when needed, on the other possibility. Large values of F seem also consistent with the MSSMJj where, for large values 



^In the case of the MSSM the Yukawa couplings enhancement is certainly not the only effect which might increase the rate. 
Since the particle spectrum changes completely from the case of the MSM new diagrams will contribute to the rate producing a 
further enhancement. In particular one should perhaps assume that the right electron number is now shared between electrons 
and selectrons by supergauge interactions. Therefore it will be necessary to consider also processes which change selectron 
number. 
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of tan/3 (giving the ratio of the two v.e.v. of the two doublets), the right electron Yukawa couplings can be enhanced 
by a factor (cos/3)^^. 

Since for T > Tc the fermion number seats both in fermions and in hypermagnetic fields (which can carry ferniionic 
number) we can expect that for T < Tc the fermion number sitting in the hypermagnetic fields should be released 
in the form real fermions because the ordinary magnetic field (surviving after the phase transition) does not carry 
fermion number. This argument holds provided the EWPT is strongly first order as it could occur in the case of the 
MSSM also for large values of tan /3 [Q . Notice that the presence of a hypermagnetic field might affect the dynamics 
of the phase transition itself |9||lq|. It was recently pointed out that, if we have a homogeneous hypermagnetic field at 
the time of the EWPT |n9| , the cross-over region m^ of the phase diagram is still present for sufficiently large values 
of the Higgs boson mass. 

Thus, for r > Th, and consequently, in the context of the MSSM, the density of Chern-Simons number present 
before the phase transition can be released in fermions which will have some chance of not being destroyed by 
sphalerons processes H and then the density of baryonic number can be related to the integrated anomaly in the 
usual way |lQ,Ol . The B AU can then be expressed directly in terms of the hypermagnetic helicity as 



UB,^ , a' UfHY-yxnYTMo 

(X,tc) = — — -;t- 



s Inac s 1 + 1 -H J c 

[s = {2/4:5)Tr'^NeffT^ is the entropy density and N^ff is the effective number of massless degrees of freedom (106.75 
in the minimal standard model)]. In the minimal standard model (MSM) T ^ T{Tr/Mo) where Tr ~ 80 TeV. If 
r > Th the hypermagnetic field has to be sufficiently weak and/or T is quite large. For example in the minimal 
supersymmetric standard model (MSSM) the right electron Yukawa coupling is larger than in the MSM case so that 
Tfj can be larger by a factor of the order of 10'^ for value of tan (3 ~ 50. If the hypermagnetic background is sufficiently 
intense, then, Th > T. This last case can arise in the MSM. If we suppose that Th > T we run, however, into a 
contradiction among our assumptions. In fact if we are just in the context of the MSM the phase transition cannot 
be strongly first order for Higgs boson masses larger than the W boson mass. This result is not crucially modified 
by the presence of a strong hypermagnetic background |19 1^. Therefore, if Th > T the BAU possibly generated 



for T > Tc hy our mechanism will be washed out by a subsequent stage of thermal equilibrium. However, even if 
we would assume that, thanks to some presently unknown phenomena, the EWPT would be strongly first order in 
the MSM, we will be able to show in the next section that the BAU obtained in the case Th > F is too small to be 
phenomenologically relevant. 

In the case where F > F^ we have that the BAU can be expressed as 



There are two (physically complementary) situations where Eqs. (6.17)-(S.18) can be applied. The first one is the 
case of a stochastic hypermagnetic background. In this case {hb/s) is strictly zero but its fluctuations (i.e. ((ns/s)^)) 
are non vanishing jll| . On the other hand in our present case the hypermagnetic field distribution is not topologically 
trivial because of the action of the pseudoscalar oscillations and, therefore, (ub/s) ^ 0. 

VII. PHENOMENOLOGICAL CONSIDERATIONS 

As a preliminary exercise let us consider the case of a Chern-Simons wave configuration Mj^ (see also p|). Suppose 
that at some moment of time the universe is populated by parity non-invariant configurations of the type discussed 



in Eq. (2.1). using the hypermagnetic field configuration obtainable from eq. (2.1) into Eq. (3.17) we get that, in the 



<^-.«=.^(?) (I) (f )«'('■)-»■"(!) (w). <-) 



limit F > Tb (wc work, for simplicity, in flat space), the BAU can be written as 

{!^\ (^ tA ^ -^^ (m fh\ fMl\ ^^(tA ^ ini" fh\ (H 

As we can see, a substantial BAU can be achieved provided the primordial value of the initial hypermagnetic field 
is large with respect to vacuum fluctuations. Suppose in fact that T-p is only built up from vacuum fluctuations. 
Then Ti^ ^ k^. Taking now into account that, at the electroweak scale, the conductivity momentum is, roughly, 
k„ ~ 10~^ Tc we have that the obtained BAU for the highest possible mode which could survive in the plasma is 
of the order of 10~^^. Too small to be relevant. In the framework of this example we see that there exist always a 
minimal hypermagnetic energy density which is compatible with the BAU. The reason for this fact is the existence of 
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a finite difFusivity scale. In our case the smallest magnetic energy density we can consider is given by Tp ^ 10~^'*T^. 
For smallest values the BAU cannot be reproduced not even taking the maximal frequency of the spectrum, namely 
fco ^ fccr- Different Chern-Simons condensates will have, of course, different numerical features and in the following we 
will try to understand if the presence of dynamical pseudoscalar can give a dynamical justification of these preliminary 
considerations. 

We will address, separately, the case of amplification occurring during inflationary phase and during radiation 
dominated phase. If the pseudoscalar oscillations take place in an inflationary phase, the production of Chern-Simons 
condensates from the vacuum fluctuations of the associated hypermagnetic fields is possible. If, on the other hand, 
the oscillations will occur in the radiation dominated epoch there is no chance of getting interesting phenomenological 
values starting from vacuum fluctuations. The reason for this behavior is very simple: if the pseudoscalar oscillations 
arise during an inflationary epoch the amplification of the hypermagnetic fields is not damped by the finite value of 
the conductivity. 
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X = log(m/H . ) 
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x= log(m/ H. ) 



FIG. 5. We illustrate the result obtained in Eq. (7.5). The two dot-dashed lines in the plots correspond to Hi ~ 10~ Mp 
and Hi ~ 10~'^^ Mp. The physical region of our parameter space is within the two dot-dashed straight lines since the maximal 
(inflationary) curvature scale cannot exceed 10~^Mp (as required by the measurements of the CMB anisotropy) and cannot 
also go below lO""^"^ Mp (where the EWPT takes place). Moreover, we would expect the maximal inflationary curvature scale 
to be larger than the curvature scale of the electroweak epoch since we want, by that time, the Universe to be dominated 
by radiation. In the left picture we report, in full lines, various curves leading to a BAU of the order of 10~^ for diflterent 
values of the coupling constant of the hypercharge field to the pseudoscalar. In the right picture we fixed c = 0.01 but we left 
z — logj^Q (tpo/M) free to change from to 2. In order to have a reasonable BAU we have to stay within the dashed lines but 
also above the thick lines. We can clearly see that this can be achieved not only in the case of large x and small y, but also 
for X ~ 1,2 and y ~ —6,-10 provided either the coupling constant is not too small and/or the initial amplitude of oscillation 
is larger than M (i.e. z > 1). In this and in the following two figures we always took N^ff = 106.75 and ctq = 70 since we 
assumed our hypermagnetic modes to re-enter in the symmetric phase of the standard electroweak theory. 



A. BAU from Hypermagnetic Knots Generated at the end of Inflation 



Eqs. (6.17)-( p.l8 ) can be explicitly evaluated using the form of our Bogoliubov coefficients given in Eqs. ( 4.13| )- 
(4.14) with the result that for co ^ i^mux the mean value of the hypermagnetic helicity is 



(Hy-V X Hy) = 



<W-)(.^) ,~n^^f 



(7.2) 



Inserting Eq. (7.2) into Eq. ( |6.18 ) and taking into account that the maximally amplified frequency toiti) ^ w,nax(^i) 
to(c/2)(V'o/-M), red-shifted at the electroweak epoch, reads 



^maxv^c 



C / V'O \ / "^ 



\/5<"'"'"'^- 



2\M J \Hi 
we can obtain the explicit form of the logarithmic spectrum of the BAU, namely 



(7.3) 



d{nB/s) ^ g'^C' AbnfiN.ff)^^ Mp /y^pX V m \ V K, 
dlogcu 512 crp7r6 Tc \M J \HiJ \Mp 



■g'=(7?T)(lf)e-2(^5r)' 



(7.4) 
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We want now to compare Eq. (7.2) with the required value of the baryon asymmetry. Defining x = log^o {'m/H., 
y = l o gin (g j/Mp) and z 
(|6l7|)-(|6Tl|), that 



logj^o ii^o/M) we have that, if we want nB/s>10 , we have to demand, from Eqs 



y > -8.5 - — logio 



N, 



eff 



logio 0-0-2 log 



10' 



2{x 



,)_l£lO(-+^)log^^e, 



(7.5) 



where ctq = o-c/T ^ 70-100 |^l|. Eq. ( |7.5| ) illustrated in Fig. ^and ||. In the thick lines we represent curves of constant 
BAU (of the order of 10~^°) as function of the two parameters of the model, namely the maximal inflationary scale 
and the mass of the pseudoscalar. We can see that for inflationary curvature scales Hi ~ 10~^-10~^*' Mp the mass has 
to be of the order of 10^ Hi. At the same time, provided the parameter space does not hit the electroweak curvature 
scale (lower dot-dashed line in Fig ^), the inflationary scale can be reduced and then, for Hi ^ 10~^^-10~^*' Mp, we 
should require masses of the order of lO'^-lO^ Hi in order to achieve the required value of the BAU. 




FIG. 6. At the left we illustrate a surface of constant BAU for a coupling constant of the order of c - 
X, y and z. At the right we plot the same quantity [i.e. Eq. (7.5)] but for c = 0.001. 



0.01 as a function of 



If, initially, the field is not in its vacuum state the hypermagnetic helicity is given by 



(iJy • V X Hy) 



H-4(Wmax)re' 



■2c(^)(^) 



(7.6) 



Inserting Eq. ( |7.6D into Eq. ( |6.18| ) the logarithmic BAU spectrum will depend explicitly on r(wniax) which is the 
critical fraction of energy density stored in the initial (topologically trivial) hypermagnetic field distribution. 

In order to be compatible with value of the BAU obtained in the context of the homogeneous and isotropic BBN 
scenario we have to require 



2/^-46. 



- 2 ^°SlO ^""ff ~ ^ ^°SlO CTO - 2 logio ^(^max) " 2 log^, c - 2{x ^ z) - 2 c 10"'+"' logio ' 
. (7.2)-(7.5) together with Eq. (7.7) we see that in order to achieve the observed 



(7.7) 



By looking at Eqs. (|7.2| )-( |7.5|) together with Eq. (|7.7| ) we see that in order to achieve the observed value of the 
BAU with reasonable values of the pseudoscalar mass in curvature units (i.e. m/Hi ~ 10^) we have to postulate an 
inflationary phase at curvature much higher than the electroweak epoch Hi ^ 10~^'^Afp. If, initially, the hypercharge 
modes are not in the vacuum we can see that if m>Hi and Hi ~ 10~^*'Mpthe observed value of the BAU can 
be achieved provided r(ci;max) ~ 10^^". This aspect is illustrated in Fig. |7[ By lowering the energy density of the 
primordial stochastic background we are pushed towards higher inflationary scales and/or higher pseudoscalar masses. 
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x= log( m/ H . ) 
FIG. 7. We illustrate Eq. ( |7.7[ ) in the case of 2 = and r(LJmax) ~ 10"^". The acceptable region of parameter space lies, 
as usual, between the two dot-dashed lines and above the full (thick) lines. If the initial state of the hypercharge fields is the 
vacuum the maximal inflationary scale compatible with a BAU of the order of 10"^" gets lower if compared with the case 
reported in Fig. BI We also illustrate (full thin lines) the same case reported in the thick lines but for a lower value of r(tJmax), 
namely, 10~^^. 



Therefore, the generation of hypermagnetic helicity is quite likely if the pseudoscalar oscillations take place during 
an inflationary phase for a wide range of masses. If we want the BAU generation to occur at very high curvature 
scales (i.e. Hi ^ lO-^-lQ-^" Mp) we are also led towards high mass scales (i.e. m>lO'^Hi ^ 10^-10^^ TeV). The 
exponential amplification appearing in the Bogoliubov coefficients is responsible of the sharp cut in the plots reported 
in Figs. I3-II. Thus, if we take the curvature scale to be sufficiently low (say, for instance, 10~^^Mp), then we can 
get a sensible BAU for masses m>10'^Hi, i.e. masses of the order of few hundred GeV at the time the amplification 
occurred. Notice that the inclusion of stimulated emission goes exactly in this direction, namely it allows, more easily, 
low curvature scales. 

We finally notice that, from Eq. (7.3), the the maximal amplified frequency (red-shifted at the electroweak scale ) 
is always smaller than the diffusivity frequency at the electroweak scale (i.e. uJa{tcw) ^ 10~^ T'cw)- This guarantees 
that the obtained BAU will not be washed out by hypermagnetic diffusion. 



B. BAU from Hypermagnetic Knots Generated During a Radiation Epoch 



We are now going to consider the phcnomcnological implications of pseudoscalar oscillations occurring in a radiation 
domin ated phase. In order to get a more explicit expression of the BAU in this case let us insert Eq. ( |3.5[ ) into Eq. 
( 3.17 ). Since Wmax/'-^o- < 1 we have that the BAU can be written as 



ub 45n/ ,A^p 



(7.8) 



It is interesting to notice that, as far as the maximal frequency is concerned, the BAU does not depend upon the 
pseudoscalar mass but only upon the pseudoscalar coupling constant. Let us estimate, very roughly, the size of the 
BAU emerging from Eq. (7.8). Suppose that N^fj = 106.75 (same particle content of the standard model) and that 
c = 0.01. Thus for ctq ^ 70 we get that m/s ~ 10~^r(a;max)- Thus the BAU can be reproduced if r(a;i„ax) ^ 10"^. 
Therefore, the presence of an initial background of hypermagnetic fields is essential as in the case of the Chern-Simons 
wave discussed at the beginning of the present Section. Suppose in fact that that A{ll}) = 1. In this case the initial 
stochastic background is formed by vacuum fluctuations. Therefore r(ajmax) = N~}A(jJraax/T)'^ ^ 10^'^^ giving a BAU 
of theorder of lO-'^o. 

The logarithmic variation of the BAU as a function of the parameters of the model can be expressed as 



y^ 



-2.4 + logio c - logio o-o + logio(-^ 



■) + 



X 



(7.9) 



where now x = logj^g ''('^max) and y = logj^p (j>'b/s). From this expression illustrated in Fig. ^ we can draw two 
conclusions. The first one is that the BAU possibly generated during the radiation dominated epoch depends upon 
the pseudoscalar coupling but not upon the pseudoscalar mass. This happens because we estimated the BAU arising 
in the case of the maximally amplified frequency whose actual value is dictated by Eq. (5.4). The smaller is the 
coupling the larger has to be the primordial hypermagnetic background. 
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FIG. 8. We illustrate the logarithinic variation of the BAU computed in Eq. (iT^) for different sets of parameters. More 
precisely we have c — 0.01, ctq — 70 (full thick line), c = 10~^, o"o = 70 (dot-dashed line), c — 10~*, cro = 100 (full thin line). 



The second observation is that the possible explanation of t he B AU through Eqs. (7.8)-( |7.9| ) depends crucially 
upon the strength of the EWPT. In fact, in order to derive Eq. (7.8) we demanded the rate of right electron chirality 
flip processes (i.e. F) to be larger than F^j. As we discussed, this hypothesis is certainly not forbidden in the case of 
the MSSM. One could also wonder if the MSM physics would be enough to implement our scenario. Let us assume. 



for a moment, to be in the framework of the MSM. Then, as argued in |lOyil| ], F// > F. The principal problem with 
this assumption is the strength of the EWPT. In fact we know that even in the presence of a strong magnetic field 
|l9| the typical cross-over behavior ||l^ of the phase transition does not change Q. 

However, even if the phase transition would be strongly first order our mechanism would be ineffective. Indeed, let 
us compute the BAU in the case where V h > F. We have that the BAU can be expressed as [M 



s 



If Hy -"^ X Hy. FMo 



2'Kac s 



Th 



r2 



(7.10) 



The expectation value appearing in the previous equation involves a highly non local operator which can be evaluated 
for scales larger than the diffusivity scale with the result that m^ 



,Hy -V X Hy. (Hy -V X Hy) 



\Hy? 
Thus, we get that the BAU can be expressed as 

s 



0.04 c 






{\Hy? 



Tr 

Mo 



(7.11) 



(7.12) 



where Tr ~ 80 TeV is the right-electron equilibration temperature |Q . For typical values of the parameters (the 
same we used in the case F > Tfj the BAU is of the order of 10~^^. We would be tempted to say that for large Tjj the 
BAU could be recovered. Unfortunately this assumption would contradict our initial hypotheses since it would mean 
that we go beyond the MSM and that, consequently, we would have Th > E. Notice that, as previously discussed P2[ , 
the results of our analysis are not in agreement with the statement that the BAU can be generated from dynamical 
pseudoscalars in the MSM, for hypermagnetic modes of the order of fcmax ~ T ^^ and independently of the specific 
value of the perturbative rate of right electron chirality flip processes. 

We want to observe, finally, that some mechanisms (like the breaking of conformal invariance) operating at some 
higher energy scale Q can generate a topologically trivial hypermagnetic background whose energy density is small 
in critical units but anyway quite large if compared to the energy density of the vacuum fluctuations. Since we 
are at high conductivity, the dissipation of the seeds is protected from being washed out by the hypermagnetic flux 
conservation. Then, at some temperature (in the TeV range) pseudoscalar oscillations will begin. At finite conductivity 
the pseudoscalar oscillations do not crucially amplify the initial hypermagnetic background but they twist its topology 
leading to a non vanishing Chern-Simons number density and, ultimately, to the BAU. 



■^More precisely, for sufficiently small hypermagnetic backgrounds (i.e. \Hy\/T^^0.3) the phase transition is stronger but 
it turns into a cross-over for sufficiently large Higgs masses. It would be interesting to understand if for sufficiently strong 
hypermagnetic backgrounds the Ambjorn-Olesen phase can be generated |45| 
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VIII. PARITY BREAKING FROM AXIONS? 

The same mechanism studied in the case of HK can be indeed exploited in order to produce ordinary magnetic 
knots for temperatures smaUer than 100 GeV. The natural candidate would then be, in this case, the ordinary QCD 
axion. We will start by recalling the possible implications of global parity breaking in ordinary MHD and we will 
then show that the amount of global parity breaking induced by axion oscillations is indeed quite minute. 

A. Global Parity Breaking in MHD 

If, at some scale ^ {H ■ V x H) ^ 0, than we also have that the Alfven (stationary) flow will be tangled in the 
same way the related magnetic field is tangled since {v ■ uj) — {v ■ W x v) j^ 0. This last condition together with the 
incompressible closure \/ ■ v — implies that also the velocity field is formed by a collection of closed loops whose 
degree of knottedness can be estimated from v-uj. Notice that in the limit of negligible thermal and magnetic dffusivity 
both the magnetic and the velocity field can be treated with topological techniques 0,0. When the thermal and 
magnetic diffusivity are adiabatically switched on (or, more physically, when we approach the magnetic and thermal 
diffusivity scales) various dissipative phenomena arise and turbulence sets in. We do not want even to attempt to 
address the very rich structure arising in MHD turbulence. However, we simply want to notice that there are crucial 
differences if the turbulence arises starting from a parity invariant or from a parity non invariant configuration of the 
bulk velocity field MH . It has been actually observed long ago that EtI , apparently, cosmic turbulence occurs usually 
in rotating bodies with density gradients. This results in a violation of the mirror symmetry of the plasma since it 
is possible, by taking the scalar product of the vorticity with the density gradient, to construct a pseudoscalar. The 
relevance of (global) parity breaking in MHD can be also appreciated by looking at the turbulent dynamo mechanism 
P,E^ . The idea is that in astrophysical and cosmological plasmas the correlation scale of t he vel ocity field is normally 



shorter than the correlation scale of the magnetic field. Under this assumption from Eq. (6.1C) it possible to obtain 
a dynamo equation by carefully averaging over the velocity field over scales and times exceeding the characteristic 
correlation scale and time tq of the velocity field. In this approximation the magnetic diffusivity equation can be 
written as 

?^ = a(V X i?) + -V'i?, a = -^{v-^xv}, (8.1) 

OT (7 O 

where a is the so-called dynamo term. In Eq. (|8.l|) H is the magnetic field averaged over times larger than tq, which is 
the typical correlation time of the velocity field). We can clearly see that the crucial requirement for the all averaging 
procedure we described is that the turbulent velocity field has to be "globally" non-mirror-symmetric. If the velocity 
field is parity invariant (i.e. no vorticity for scales comparable with the correlation length of the magnetic field), then 
the dynamics of the infrared modes is decoupled from the velocity field since, over those scales, a = 0. Therefore, for 
our purposes the a term will be a good measure of the degree of knottedness associated with the velocity field. 

Needless to say that that there is a puzzle associated with the turbulent dynamo mechanism [^j. On one hand 
we see that our galaxy has an a term which can be roughly estimated to be ac ^ 0.8 x 10^ cm/sec in the interval 
3 Kpc<r<12 Kpc (it is assumed that \uj\ ~ 10~^^ sec~^) H. On the other hand we have to postulate that some alpha 
term had to be present also at the time of the protogalactic collapse (occurred around tgc ~ 10^^ sec). This is because 
we want the dynamo term to amplify, by differential rotation the primordial magnetic seeds up to the observed value 
of themagnetic field observed in spiral galaxies. Thus, the parity breaking of the galactic turbulence seems to be more 
a consequence of some primordial (global) parity breaking than a theory of the galactic dynamo action. 

If one wants to postulate some global parity breaking in the plasma the obvious question which arises is when this 
phenomenon could occur in the early Universe. It should certainly be useful for the dynamo action to have some sort 
of global parity breaking at the decoupling epoch {tdec ^ lO^^sec) which is not so far (for cosmological standards) 
from the time of the protogalactic collapse. In principle, the parity breaking could occur even before but the only 
problem could be, in this case, that many of the known phenomena occurring in the early Universe are normally 
discussed in terms of a mirror-symmetric plasma. For example the present BBN calculations do not include, to the 



best of our knowledge, possible effects of global parity violation |48 49] and we would guess that to have global parity 
breaking might be a problem for nucleosynthesis. 



We denote with H the ordinary magnetic field as opposed to the hypermagnetic field. 
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The axion oscillations starting at a temperature of the order of 1 GeV can certainly produce a source of global 
parity breaking. The question is how large is the induced a term in the primordial plasma. Our aim is now to estimate 
the order of magnitude of the parity breaking induced through primordial electromagnetic helicity. 

B. Parity Breaking induced by ordinary Magnetic Knots 

As it is well known in axion models |^,0 an extra (global) Upq{1) symmetry complements the MSM. This 
symmetry is broken at the Peccei-Quinn scale Fa and leads to a dynamical pseudo Goldstone boson (the axion) which 
acquires a small mass because of soft instanton effects at the QCD phase transition. If an axionic density is present in 
the early Universe , bounds can be obtained for the Peccei-Quinn symmetry breaking scale. These bounds together 
with astrophysical bounds (coming mostly from red-giants cooling and from SN 1987a) leave a window opportunity 
10i*^GeV <Fa< IQi^GeV Q. 

There are different realizations of axion models. In some models ||2J,Q quarks and leptons do not carry the Peccei- 
Quinn charge, whereas in other realizations the usual standard model leptons are not singlet under the Peccei-Quinn 
symmetry pJpS]. We will not enter into the details of these two models since what matters for our consideration is, 
at least in the very first approximation, how the axion couples to the anomaly. 

There are also different cosmological realizations of axionic models. If an inflationary epoch takes place prior to 
the usual radiation dominated regime, the Peccei-Quinn phase can be aligned by today |54| . In the absence of an 
inflationary phase it is likely that as the Universe passes through the Peccei-Quinn phase transition a network of 
cosmic strings will be generated via the Kibble mechanism ||55| . We will assume that an inflationary phase took place 
weh before (H < lO'^Mp) the onset of the EWPT. 

If the temperature of the Universe exceeds the critical temperature of the EWPT the effective action describing 
the interaction of axions with the Standard model gauge fields can be written as 



5, 



eff 



I (ai^-^'^" + <^^w£,K.^^'' + ^^- 3&^-^^'' 



■2) 



where (7^, g and g' are, respectively, the SUc{S), SUl{'2.) and C/y(l) coupling constants [c^w ^nd Ctpy are numerical 
constants of order 1 which depend upon the specific axion model). 
For T < Tr the effective action becomes 



Seff = I d^ 



i' f 9. 



2 






KV 32712^"^" ■ """'32^' 



(8.3) 



where F^j^^, is the ordinary (electromagnetic) gauge field strength. 

There are important differences between the axion evolution for temperatures larger and smaller than the EWPT 
temperature. In fact |56| , for T > Tc the axion evolution equation is not only modified by the Hubble damping factor, 
but also by damping induced by the QCD sphalerons which are not suppressed at high temperatures. This effect 
leads to an equation for tp which differs from the one introduced in the previous section and which can be written as 

i) + {3H + -f)ij = (8.4) 

where we dropped for simplicity the potential term. It turns out that |pq] 

^ = ;^^ ^ -^ (8-5) 

(where as — gl/An). Thus, if Fa > lO^GeV we have to conclude that the sphaleron induced damping dominates 
over the damping produced by the expansion of the Universe (i.e. ~ T'^/Mp). The axion time evolution will then be 
(approximately) suppressed as "0 ex e~^* ~ e~'^'^ in a radiation dominated Universe (recall that t ~ t^). Therefore, 



we have to accept that in Eq. (2.18) the polarization mixing will be exponentially suppressed as time goes by. The 
sphaleron induced damping does not exclude other possible consequences associated with the axions at the EWPT 
like the possible existence of strong CP violating effects [M 

Notice that this argument does to not apply to temperatures T <^Tc when the plasma has already cooled sufficiently 
for 7 to have turned off. Therefore, for our purposes it seems certainly theoretically more plausible to consider the 
possible occurrence of coherent axion oscillations at temperatures of the oder of 1 GeV when the electroweak symmetry 
is broken. 
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In can be, in principle, dangerous to have production of parity non invariant configurations of the ordinary magnetic 
field right before the onset of BBN. Suppose that the axion oscillations start at T ~ 1 GeV with a typical mass of 
j^Q-ig QgY_ Therefore we have to consider the ordinary axion oscillations coupled to the evolution of the ordinary 
vector potential at finite conductivity 






ctA'± 



fc^Tfc- 



27r Fn 



A± = 0. (8.6) 



At high temperatures T 3> ^qcd the axion is massless, but, at lower temperatures it develops a mass due to QCD 
instanton effects. The effects of instantons is highly suppressed at high temperatures, however, the temperature 
dependence of the axion mass is approximately given by [ p8| 

m(T)^ 0.1 (^^) mo, (8.7) 

where toq is the axion mass at zero temperature. 

Axion oscillations between a temperature of the order of 1 GeV and T ^ 10 MeV (i. e. around the nucleosynthesis 
epoch) can generate electromagnetic helicity. Suppose, for sake of simplicity that the axion oscillations will start with 



a typical mass f the order of 10 ^^ GeV. Then, according to Eq. (8.6), we have that the maximal amplified frequency 
and the diffusivity frequencies are of the order of 



^^^, ..^.K^,/^T. (8.8) 

47r Fa Mo' V aem \l Mp "• ' 



where T is of the order of 1 GeV. We can solve Eq. (B.E) similarly to what we did in hypermagnetic case. In the case 



(^max < ^CT the magnetic helicity produced via axion oscillations will be of the order of 

10-22 (8.9) 



{H-WxH) ^ aemAV' T, ^^_^^ 



apc 47r Fa Mq 



(notice that (Tq ^ a^^^^ and Ta ^ 1 GeV). Eq. (S.9) provides then a rather small parity breaking and certainly irrelevant 
for the dynamo mechanism. 

IX. CONCLUDING REMARKS 

In this paper we investigated the production of hypermagnetic knots. We found that HK can be efficiently produced 
provided some topologically trivial configuration of the hypermagnetic field is already present for scales larger than the 
hypermagnetic diffusivity scale. The BAU can be successfully reproduced, but not in the MSM. The reason is twofold. 
On one hand the phase transition is too weak (even taking into account the possible presence of the hypermagnetic 
background). On the other hand the (perturbative) rate of the slowest processes in the plasma is too small. In the 
MSSM the BAU generation seems more probable and our mechanism can be viable, provided the typical scale of the 
HK exceeds the hypermagnetic diffusivity scale prior to the onset of the EWPT. Inspired by the axial Higgs we gave 
examples of specific scenarios where this mechanism can be implemented. We also discussed the possible production 
of magnetic helicity by ordinary QCD axions: we found that its value is too small to be phenomenologically relevant. 
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APPENDIX A: EXAMPLES OF HYPERMAGNETIC KNOTS FROM TOPOLOGICALLY NON TRIVIAL 

CONFIGURATIONS OF MAGNETIC FIELDS 



In Section II we presented some topologically non trivial configurations of the hypercharge field [g||9|] whose main 
feature is to have a Chern-Simons number density which is time dependent and uniformly distributed in space. It 
is also possible to construct hypermagnetic knot configurations with finite energy and helicity which are localized in 
space and within typical distance scale L^ . Let us consider in fact the following configuration in spherical coordinates 

2Bo COS0 



ye{n,t 



y^{n,6)^- 



^Ls [7^2 + 1] 

2Bo sin 6* 
2Bo nU sm0 



2 ' 



Trie 



(B.l) 



[n^ + iy 

where TZ = r/Ls is the rescaled radius and Bq is some dimensionless amplitude and n is just an integer number whose 
physical interpretation will become clear in a moment. The hypermagnetic field can be easily computed from the 
previous expression and it is 

4i?o n cos 9 



nr{n,e) = - 



Hein, 



n4,in,i 
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[7^2 + l]- 
7^2-l 



7rL2 
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TT LI 



[R? + l\ 



The poloidal and toroidal components of Ti can be usefully expressed as Tip — TirSr 
Chern-Simons number is finite and it is given by 



Ncs = 



9' 
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(B.2) 
Tiees and Tit = Ti^e^. The 

(B.3) 



In Eq. (B^) the integration is perfectly convergent over the whole space and, since the field is Tiy = in any part of 
dV at infinity, this quantity is also gauge invariant. Thus, we can see that the integer index n labels the number of 
knots of the configuration. We can also compute the total helicity of the configuration namely 

256 Bl n [°" TZ^dU hBln 

" ' - ^ (B.4) 



Wy • V X T-iyd-^x = 
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L? 



There is an important difference between Eqs. (B.3) and (B.4). In Eq. (B.3) the integrand is not g auge invariant (but 
the integrated Chern-Simons number density is gauge-invariant), on the other hand in Eq. (B.4) also the integrand 
is gauge- invariant. We can compute also the total energy of the field 



E 



d^x\HY? 



Bl 
2Ls 



{n^ 



1) 



(B.5) 



and we discover that it is proportional to n^ . This means that one way of increasing the total energy of the field is 
to increase the number of knots and twists in the flux lines. We can also have some real space pictures of the core of 
the knot (i.e. TZ = r/Lg < !)■ In Fig. g wc plot the cartesian components of the hypermagnetic field 



ny{x) = - 



4So 



2X 



y^ + z^\ 

V2nyZ 



HAX) 



^Li \i^x'^ + y^ + z'^Y 

4Bon[X^+y^-Z^-l] 
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[1 + ^-2+^^2+22]^ 



(B.l 



26 



in the case of n = and n = 5 in terms of the rescaled coordinates X = x/Lg. y = y/Lg, Z — z/Lg 




FIG. 9. We plot the hypermagnetic knot configuration in Cartesian components [Eqs. (B.6)] for the cases n — (zero 
helicity) and n = 5. The direction of the arrow at each point represents the tangent to the flux fines of the magnetic field, 
whereas the length of the vector is proportional to the field intensity. We notice that for large n the field is concentrated, in 
practice, in the core of the knot (i.e. r/Ls < 1). 



APPENDIX B: EXPECTATION VALUE OF THE HYPERMAGNETIC HELICITY 



The aim of this Appendix is to compute the expectation value of the hypermagnetic hehcity which represents a fair 
gauge-invariant estimate of parity breaking induced by the presence of dynamical psudoscalars in the hypermagnetic 
background. 



Let us expand the hypercharge field operator according to Eq. (3.15) but with slightly different notations which 
turn out to be useful for the calculation at hand: 



Yout {x, t) 



V2^J 7(2^ 



Yp.out{k,T)efie'''-^ + rJ_„„,(fc,T)e^e 



-ik-x 
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(where, as usual, /3 = +, — ). 

Therefore the expectation value of the hypermagnetic helicity will be 



(0|J?y -Vxi/ylO) 



1 



(fk 



dPp 



fV 



(y+(fc, r)ft (p, r))e'^^-^-^ - {Y_{k, t)y1{p, r))e'(>^-P^-^ 



2 J y/{2^J 7(2^3 

+ (r;(fc,T)y+(p,T))e-^(^~p>"^- (r!(fc,T)r_(p,r))e-'(*^"p^-"^ 



(B.2) 



(we dropped the "out" subscript for the mode functions). The various correlators appearing in the previous expression 
can be evaluated, after some algebra, using the explicit form of the Bogoliubov transformation given in section II: 
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From these last expressions we get Eq. (5.16). With the same technique other expectation values of parity non- 
invariant operators can be computed like 
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Notice that this second gauge-invariant operator can be also used in order to measure the amount of parity breaking. 
However, in the plasma it seems that the first operator computed in this Appendix is more suitable. Indeed, at finite 
conductivity, inhomogeneous hypermagnetic fields lead to the presence of induced (Ohm) hyperelectric fields of the 
order of Ey ~ cr~^V x Hy 
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